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From the Aharonov-Bohm effect to general relativity, geometry plays a central role in modern 
physics. In quantum mechanics many physical processes depend on the Berry curvature. However, 
recent advances in quantum information theory have highlighted the role of its symmetric counter- 
part, the quantum metric tensor. In this paper, we perform a detailed analysis of the ground state 
Riemannian geometry induced by the metric tensor, using the quantum XY chain in a transverse 
field as our primary example. We focus on a particular geometric invariant - the Gaussian curvature 
- and show how both integrals of the curvature within a given phase and singularities of the curva- 
ture near phase transitions are protected by critical scaling theory. For cases where the curvature is 
integrable, we show that the integrated curvature provides a new geometric invariant, which like the 
Chern number characterizes individual phases of matter. For cases where the curvature is singular, 
we classify three types - integrable, conical, and curvature singularities - and detail situations where 
each type of singularity should arise. Finally, to connect this abstract geometry to experiment, we 
discuss three different methods for measuring the metric tensor, namely via integrating a properly 
weighted noise spectral function and by using leading order responses of the work distribution to 
ramps and quenches in quantum many-body systems. 



Understanding the geometry and topology of quantum 
ground state manifolds is a key component of modern 
many-body physics. For example, it has become stan- 
dard practice to characterize topological phases by their 
Chern number, defined as the integral of Berry curva- 
ture over a closed manifold in parameter sp ace. Exam- 
ples of this include the quantum Hall effect!^, topological 
insulators^', integer and half integer spin chains^!, and 
many others. Non-zero Berry curvature is typically as- 
sociated with broken time reversal symmetry, either ex- 
plicitly by external coupling to a time reversal breaking 
field or implicitly by splitting the ground state manifold 
into diffe rent sectors, each of which breaks time reversal 
symmetry^! 

However, in addition to the Berry curvature, which de- 
scribes the flux of Berry phase within the ground state 
manifold, another important quantity is its symmetric 
counterpart - the quantum (Fubini-Study) metric tensor 

- which describes the absolute value of the overlap am- 
plitude between neighboring ground states^. This metric 
plays an important role in understanding the physics of 
quantum many-body ground states^^and is at the heart 
of current research in quantum information theorjW^. 
For instance, the diagonal components of the quantum 
metric tensor are none other than fidelity susceptibilities, 
whose scaling in the vicinity of quantum phase transitions 

- including topological phase transitions - is an object of 
great interesB^ES. 

The purpose of this paper is to understand the quan- 
tum geometry of a simple model, the spin- 1/2 XY chain 
in a transverse field. For this integrable model, we solve 
the geometry and topology of the ground state manifold 
as a function of three parameters: transverse magnetic 
field, interaction anisotropy, and spin rotation about the 
transverse axis. Using a standard trick from Ricmann 
geometry, we analyze the three-dimensional metric by 



taking a series of two-dimensional cuts. For cuts along 
which the Riemannian manifold is regular, we identify 
the shape of the manifold and show that the shape of 
each phase is protected against symmetry-respecting per- 
turbations by critical scaling theory of the metric tensor. 
For cuts along which the manifold is singular, we identify 
and classify the singularities. As with the other cuts, we 
demonstrate that the singularities are robust against a 
variety of modifications to the low-energy theory. We see 
three types of geometric singularities: integrable, conical, 
and curvature singularities. Finally, we detail general cir- 
cumstances where each type of singularity will arise. 

Given the importance of the quantum metric tensor in 
understanding the properties of ground state manifolds, 
it is surprising that there have been no direct experimen- 
tal measurements of the metric tensor to date. There- 
fore, at the end of the paper we discuss several different 
proposals for experimentally measuring the components 
of the metric tensor. The first method is based on the 
direct representation of the metric tensor through the 
noise spectral function, generalizing a recent proposal 
by Neupert et al.^ for measuring the metric tensor in 
non-interacting Bloch bands. The second method re- 
lates the metric tensor to a measurement of the lead- 
ing non-adiabatic contribution to the excess heat for 
square root ramps (see also Refs. flBl and [TT|) . The third 
method similar identifies the metric with leading non- 
adiabatic corrections to energy fluctuations for generic 
linear ramps. Finally, the fourth method is based on an- 
alyzing the probability of doing zero work in single or 
double quenches, which is related to the time average 
of the well-known Loschmidt echcSS. Using these tech- 
niques the full many-body metric tensor is - at least in 
principle - experimentally accessible. We note that these 
measurement proposals do not rely on many of the geo- 
metric notions discussed elsewhere in the paper, so those 
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primarily interested in measuring the metric can skip di- with the ground state manifold: 



rectly Sec. IV We also note that the metric tensor can be 
readily extracted numerically as a non-adiabatic response 
of physical observables to imaginary time rampi^, by di- 
rectly evaluating overlaps of the ground state wave func- 
tions at slightly different couplings 19 , or through numer- 
ical integration of im aginary time noise spectra of the 
generalized forcedM2U. 

The rest of this paper proceeds as follows. In Sec.[IJ we 
introduce the metric tensor and show its relation to the 
Berry connection operators. In Sec. |I A| we describe how 
the Euler characteristic, curvature, and other geometric 
invariants are obtained from this metric. As a useful ex- 
ample, in Sec. |I C| we explicitly solve for these quantities 
for the case of the integrable quantum XY model in a 
transverse field. In Sec. |ID| we show how to visualize 
the metric manifold by mapping to an isometric surface 
embedded in three dimensions. These shapes motivate us 
to define invariant integrals consisting of the contribution 
to the Euler characteristic within a given phase. We solve 
this exactly for the XY model then, in Sec. [TTJ we argue 
based on critical scaling of the metric tensor that the ge- 
ometric integrals remain unchanged for all models in the 
same universality class. To further understand the geom- 
etry of these Ricmann manifolds, we classify three types 
of singularity in the Gaussian curvature that can occur in 
the vicinity of phase transitions: integrable (Sec. Ill A I, 



conical (Sec. IIIB), and curvature (Sec. IIICl singular- 
ities. Abstracting away from the XY model, we detail 
situations under which each type of singularity should 
arise. Finally, in Sec. [TV] we discuss different methods for 
measuring the quantum metric in terms of a more tra- 
ditional condensed matter measurement of noise correla- 
tions, as well as through real-time ramps and quenches 
of the system parameters as is more relevant to isolated 
cold atom experiments. 



9nu = %X^\ 



(3) 



However, in another seminal worlPI Berry introduced the 
notion of geometric phase (a.k.a. Berry phase) and the 
related Berry curvature, which is given by the imaginary 
(antisymmetric) part of the geometric tensor: 

Ffa, = -23?[x Mt ,] = i(x?v - X^) = dpA v - d u Af, , (4) 

where — i^ol^ul^o) 1S ^ ne Berry connection within 
the ground state manifold. The Berry phase is just a line 
integral of the Berry connection or - by Stokes theorem 
- a surface integral of the Berry curvature: 



A ■ dX 



<)S 



(5) 



where dS^ v is a directed surface element. 

In this work, we will be primarily interested in the met- 
ric tensor, g^ v . One simple physical interpretation of the 
metric tensor is that it sets natural units, allowing one 
to compare different physical parameters. For example, 
if we consider the ground state manifold as a function of 
magnetic field and pressure, one can ask how one Tesla 
compares to one Pascal. In the absence of a simple single 
particle coupling, the method for scaling these quantities 
to compare them is not obvious. However, the metric 
tensor provides a natural answer by allowing one to com- 
pare the effects of these couplings on the ground state 
fidelity. Rescaling the units by the corresponding diago- 
nal components of the metric tensor - a.k.a. the fidelity 
susceptibilities - one sets natural units for different cou- 
plings. Therefore, one Tesla can be compared to one 
Pascal by comparing the " dimcnsionless" couplings after 
rescaling dX^ 



GEOMETRY OF THE GROUND STATE 
MANIFOLD 



Geometric invariants of the metric tensor 



Consider a manifold of Hamiltonians described by 
some coupling parameters A. A natural measure of the 
distance between the ground state wave functions |^>o) 
separated by infinitesimal dX is^ 

ds 2 = 1 - |(^ (A)|^o(A + dA))| 2 = Y,X»vd\»d\ v , (1) 



where % M „ is the geometric tensor: 

Xp.u = (MK d v\i>o) ~ (^o|c^|V>o)(V>o|<9i>|V ; o) , (2) 

with <9 M = gfjj-. As noted by Provost and Vallee 8 , 
this tensor is invariant under arbitrary A-dependent U(l) 
gauge transformation of the ground state wave functions. 

Strictly speaking, Eq.[T]utilizes only the real symmetric 
part of Xfivi which defines the metric tensor associated 



Underlying the classification of most topological phases 
is the fact that the Berry curvature satisfies the Chern 
theorem 22 , which states that the integral of the Berry 
curvature over a closed two-dimensional manifold M. in 
the parameter space is 2tt times an integer n, known as 
the Chern number: 



M 



F^vdS^ = 27m . 



(6) 



Physically this theorem reflects the single valuedness of 
the wave function during adiabatic evolution: imagine 
splitting M. into "upper" and "lower" surfaces. To main- 
tain single valuedness, the Berry phases obtained by in- 
tegrating the Berry curvature over the upper and lower 
surfaces can only be different by a multiple of 2tt. 

While the Berry curvature and its associated geome- 
try are certainly of great interest in modern condensed 
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matter physics, the metric tensor g^ v also plays an im- 
portant role. This tensor defines a Riemannian manifold 
associated with the ground states, and it is interesting to 
similarly inquire about its geometry and topology. In 
particular, the shape of the Riemannian manifold de- 
fines a different topological number, given by applying 
the Gauss-Bonnet theorerrPH to the quantum metric ten- 
sor: 



1 

2^ 



KdS- 



M 



dM 



kgdl 



= X(M) 



(7) 



where x(A / I) is the integer Euler characteristic describing 
the topology of the manifold Ai with metric g^ v . The two 
terms on the left side of Eq. [7] are the bulk and boundary 
contributions to the Euler characteristic of the manifold. 
We refer to the first term, 



Xbulk(-M) = r- / KdS, 



(8) 



and the second term, 



Thus, for example, the bulk Euler integral (Eq. [8]) can be 
used for classification of geometric properties of different 
phases. This geometric invariant is in general different 
from the Chern number, and can be non-trivial even in 
the absence of time reversal symmetry breaking. 

While the Gauss-Bonnet theorem has a higher dimen- 
sional generalization known as the Chern-Gauss-Bonnet 
theorem- 5 , in this work we will focus only on the two- 
dimensional version. We emphasize that the dimension- 
ality here is that of parameter space; the physical dimen- 
sionality of the system can be arbitrary. The choice of 
parameters is also arbitrary, and is usually dictated ei- 
ther by experimental accessibility, symmetry properties 
of the system, or other related considerations. Choosing 
appropriate parameters for studying geometric properties 
of the phases is therefore similar to choosing parameters 
defining the phase diagram. 



B. Defining the geometric tensor through gauge 
potentials 



Xboundary 

(M) 



1 

2^T 



dM 



kgdl 



(9) 



as the bulk and boundary Euler integrals, respectively. 
These terms, along with their constituents - the Gaus- 
sian curvature (K), the geodesic curvature (k g ), the area 
element (dS), and the line element (dl) - are geometric 
invariants, meaning that they remain unmodified under 
any change of variables. More explicitly, if the metric is 
written in first fundamental form as 



ds 2 = EdX\ + 2Fd\ 1 d\ 2 + GdXj , 
then these invariants are given bjl^H 



K = 



(10) 
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dS = y/gdX\dX2 
dl = VGdX 2 , 



(11) 



where k g and dl are given for a curve of constant Ai. The 
metric determinant g and Christoffel symbols are 



g = EG 
1 



F 



ij ~ 



g km (djg im + d l g ]m - d m gij) 



(12) 
(13) 



where g lJ is the inverse of the metric tensor gij . 

As we will see, in general the bulk and the boundary 
terms are not individually protected against perturba- 
tions for an arbitrary manifold M.. A major purpose of 
the current work is to demonstrate that if the parame- 
ter space manifold terminates at a phase boundary, how- 
ever, then not only is the sum in Eq. [7] protected against 
various perturbations, but so is each term individually. 



It can be convenient to express the geometric tensor 
through the Berry connection operators A^ = id^ asso- 
ciated with the couplings A M , which one can think of as 
gauge potentials in parameter space. These gauge oper- 
ators are formally defined through the matrix elements 



A mn 



i{rn\d,j\n) , 



(14) 



which implicitly depend on the U(l) phase choice for 
each energy eigenstate |n) at each A. If the basis de- 
pendence of A is expressed through a unitary rotation of 
some parameter-independent basis, 



|n(A)> = U nm (X)\m) Q , 
then the gauge potentials can be written as 



(15) 



The operator generates infinitesimal translations 
of the basis vectors within the parameter space. For in- 
stance, if spatial coordinates play the role of parameters, 
then the corresponding gauge potential is the momen- 
tum operator. If the parameters characterize rotational 
angles, the gauge potential is the angular momentum op- 
erator. 

The ground state expectation value of the gauge oper- 
ator is by definition the Berry connection 



Ap = W>olA*|V>o) 



(16) 



The geometric tensor is the expectation value of their 
covariance matrix: 

Xpu = {ip^A^A^ipojc 

= W[A^Au\^q) - (^ |^|^o)^o|X|Vo)-(17) 

More explicitly, the components of the metric tensor and 
the Berry curvature are expressed through the connected 
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expectation value of the anti-commutator and the com- 
mutator of the gauge potentials, respectively: 



1 



= -(ipolAfj.Au + A v A^\ip ) c 
Ffju = i(ijjv\AfjA u - AvA^tpo). 



(18) 
(19) 



Although we will be interested only in the ground state 
manifold for the remainder of this paper, we briefly com- 
ment that the definitions above can be extended to arbi- 
trary stationary or non-stationary density matrices. For 
example, with a finite temperature equilibrium ensemble, 
one can define = iTr(p t herniai[-A^) A v ]). 



In Sec. |I C| we will explicitly calculate the gauge poten- 
tials of the quantum XY chain. Here we comment on a 
few of their general properties. First, we note that gauge 
potentials are Hermitian operators. This follows from dif- 
ferentiating the identity (n(X)\m(X)} = 5 nm with respect 
to A M , or more directly from their definition in terms of 



unitaries (Eq. 15). Also, the gauge potentials satisfy re- 
quirements of locality. In particular, if the Hamiltonian 
can be written as the sum of local terms % = JY hi and A 
represent global couplings within this Hamiltonian, then 
the geometric tensor is extensive^. In particular, this 
implies that fluctuations of the gauge potentials are also 
extensive, which is a general property of local extensive 
operators. Similarly, if A represent local (in space) per- 
turbations, then the geometric tensor is generally system 
size independent, so that _A M is again local. As usual, 
various singularities - including those breaking locality 
- can develop in gauge potentials near phase transitions. 
Finally, we point out that if A M is a symmetry of the 
Hamiltonian, i.e., if the Hamiltonian is invariant under 



A„ 



Xfj + SXfj, then all gauge potentials A v are also 



invariant under this symmetry. 



C. Metric tensor of the quantum XY chain 

As our primary example, we consider a quantum XY 
chain described by the Hamiltonian 



(20) 



where J x ^ v are exchange couplings, h is a transverse field, 
and the spins are represented as Pauli matrices s x,v ' z . It 
is convenient to re-parameterize the model in terms of 
new couplings J and 7 as 



7 J <^)' 



(21) 



where J is the energy scale of the exchange interaction 
and 7 is its anisotropy. We add an additional tuning pa- 
rameter </>, corresponding to simultaneous rotation of all 
the spins about the z-axis by angle (p/2. While rotating 
the angle (p has no effect on the spectrum of H, it does 
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FIG. 1: Ground state phase diagram of the XY Hamilto- 
nian (Eq. 22 1 for = 0. The rotation parameter <p modifies 
the Ising ferromagnetic directions, otherwise maintaining all 
features of the phase diagram. As a function of transverse 
field h and anisotropy 7, the ground state undergoes con- 
tinuous Ising-like phase transitions between paramagnet and 
ferromagnet at h — ±1 and anisotropic transitions between 
ferromagnets aligned along X and Y directions (X/Y-FM) at 
7 = 0. These two types of phase transition meet at multi- 
critical points, which are described in detail in Ref. 1261 



modify the ground state wave function. To fix the overall 
energy scale, we set J = 1. 

The Hamiltonian described above can be written as 

H(h, 1 ,<p) = -j2[4 s i+i+ h - c l 

-7^^44+1+^-]-^ 4- ( 22 ) 

3 3 

Since the Hamiltonian is invariant under the mapping 
7 — > —7, cp — > tp + 7r, we generally restrict ourselves 
to 7 > 0, although we occasionally plot the superfluous 
7 < region when convenient. This model has a rich 
phase diagram ^ 6 ! 27 *, as shown in Fig. [l] There is a phase 
transition between paramagnet and Ising ferromagnet at 
\h\ = 1 and 7^0. There is an additional critical line 
at the isotropic point (7 = 0) for \h\ < 1. The two 
transitions meet at multi-critical points when 7 = and 
\h\ = 1. Another notable line is 7 = 1, which corre- 
sponds to the transverse-field Ising (TFI) chain. Finally 
let us note that there are two other special lines 7 = 
and \h\ > 1 where the ground state is fully polarized 
along the magnetic field and thus ^.-independent. Thus 
this line is characterized by vanishing susceptibilities in- 
cluding vanishing metric along the /i-direction. As we 
discuss in Sec. |IIID| such state is fully protected by the 
rotational symmetry of the model and can be terminated 
only at the critical (gapless) point. The phase diagram 
is invariant under changes of the rotation angle (p. 

Rewriting the spin Hamiltonian in terms of free 
fermions via a Jordan- Wigner transformation, % can 
be mapped to an effective non-interacting spin one-half 
modePS with 



(23) 



U k = - 



h — cos(k) 
7sin(fc)e- i<?i 



7 sin(fc)e^ 
-\h— cos(fc)] 



5 



This mapping yields a unique ground state throughout 
the phase diagram by working in a particular fermion 
parity sector^ none of the conclusions below will change 
if the other sector is chosen in cases when the ground 
state is degenerate. A more gener al analysis involving 
the non-Abelian metric tensor^Eni j s outside the scope 
of this work. 

The ground state of Hk is a Bloch vector with az- 
imuthal angle and polar angle 



tan 



7sin(fc) 
h — cos(fc) 



(24) 



To derive the components of the metric tensor, we start 
by considering the gauge operators introduced in Sec.|IB| 



Afj_ = idfj_. If we consider the transverse field h, we see 



that 
9h\gs k ) 



d h 9 k 



-sin(f) 
cos(^) 



e -#/2 



dht 



esfc) 



(25) 

The same derivation applies to the anisotropy 7, since 
changing either 7 or h only modifies 6 k and not (f>. Thus 
we find 



(26) 




FIG. 2: Equivalent graphical representation of the phase dia- 
gram of the transverse field Ising model (7 = 1) in the h — 4> 
plane (Eq. 341. The ordered ferromagnetic phase maps to 
a cylinder of constant radius. The disordered paramagnetic 
phases h > 1 and h < — 1 map to the two hemispherical 
caps. The inset shows how the cylindrical coordinates z and 
r depend on the transverse field h. 



The expressions for the metric tensor can be evaluated 
in the thermodynamic limit, where the summation be- 
comes integration over momentum space. It is convenient 
to divide all components of the metric tensor by the sys- 
tem size and deal with intensive quantities g^ v — > g^/L. 
Then one calculates these integrals to find that 



where A = {h, 7} and T^ y,z are Pauli matrices that 
act in the instantaneous ground/excited state basis, i.e., 
r fe|g s fc) = |gs fc >, Tfclesfe) = -|es fc ). Similarly, for the pa- 
rameter <f>, we find that 

A * = -\Y, [ C0S (^) T I + sin (#fcK] . 



(27) 



In terms of these Berry connections, the metric tensor 
and Berry curvature can be written as (see Sec. IB) 



9/u> — 2 (vA") -4f})c 



(28) 



In the case of the XY model, the metric tensor reduces 
to 




H 
M+i' 
\h\ 



-7 2 ^^-1+7 



- 1 



\h\<\ 
\h\>l 



l7l(l-£ 2 )' 
(h 2 -l)(h 2 -l+7 2 )3/ 2 ' 

I7l(l + l7l) 2 
\h\ 



\h\<\ 
\h\>\ 

1 



\h\l 2 

(l-7 2 )(/( 2 -l+7 2 ) 3 / 2 




0, 

-\h\l 
/i(7i 2 -l+ 7 2 ) 3 / 2 



\h\<l 
, \h\>l 



(31) 



9h h 



9hj 



'd<h 
4 V V dh 

4^ dh 5 7 



9jf 



4 ^ 

fc 



c?7 



ghcj> = g-y4> 







(29) 



Although we will not be interested in the Berry curvature 
in this work, we show the corresponding expressions for 
completeness 



k 



F-i* - 2 



It, 



A- 



86, 



h), F hy =0. 



(30) 



D. Visualizing the ground state manifold 

Using the metric tensor we can visualize the ground 
state manifold by building an equivalent (i.e., isometric) 
surface and plotting its shape. It is convenient to focus on 
a two-dimensional manifold by fixing one of the parame- 
ters. We then represent the two-dimensional manifold as 
an equivalent three-dimensional surface. To start, let's 
fix the anisotropy parameter 7 and consider the h — <p 
manifold. Since the metric tensor has cylindrical sym- 
metry, so does the equivalent surface. Parameterizing 
our shape in cylindrical coordinates and requiring that 



dz + dr 



2 , r 2 £ 



9hhdh 2 + g^dcj) 2 , (32) 
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we see that 



ih) 



dhi i 



(33) 

Using Eq. |31[ we explicitly find the shape representing 
the XY chain. In the Ising limit (7 = 1), we get 



r(h) = \ 

z{h) = ~ W 

r{h) 



z(h) = 



1 

i\h\ 

_ 7T \H\ 

8 h 



ih 



\h\ < 1, 

\h\>l 



(34) 



The phase diagram is thus represented by a cylinder 
of radius 1/4 corresponding to the ferromagnetic phase 
capped by the two hemispheres representing the param- 
agnetic phase, as shown in Fig. [2j It is easy to check 
that the shape of each phase does not depend on the 
anisotropy parameter 7, which simply changes the as- 
pect ratio of the cylinder. By an elementary integration 
of the Gaussian curvature, the phases have bulk Eulcr 
integral for the ferromagnetic cylinder and 1 for each 
paramagnetic hemisphere. These numbers add up to 2 as 
required, since the full phase diagram is homeomorphic 
to a sphere. From Fig. [2j it is also clear that the phase 
boundaries at h = ±1 are geodesies, meaning that the 
geodesic curvature (and thus the boundary contribution 
Xboundary) is zero for a contour along the phase bound- 
ary. As we will soon see, this boundary integral protects 
the value of the bulk integral and vice versa. 

We can also reconstruct an equivalent shape in the 7— (f> 
plane. In this case we expect to see a qualitative differ- 
ence for \h\ > 1 and \h\ < 1 because in the latter case 
there is an anisotropic phase transition at the isotropic 
point 7 = 0, while in the former case there is none. These 
two shapes are shown in Fig. [3] The anisotropic phase 
transition is manifest in the conical singularity develop- 
ing at 7 — 0. 

The singularity at 7 = yields a non-trivial bulk Euler 
integral for the anisotropic phase transition. To see this, 
consider the bulk integral 



dep / d 7 Vs(7,<W7,0) • (35) 



Xbuik(e) = lim 

L— Voc 



In the limit e — > + , this integral has a discontinuity as 
a function of h at the phase transition, as seen in Fig. [3] 
Thus, Xbuik = Xbuik (e = + ) can be used as a geometric 
characteristic of the anisotropic phase transition. As we 
will show in the next section, its value is Xbuik = l/\/2 in 
the ferromagnetic phase and Xbuik = 1 in the paramag- 
netic phase. This non-integer geometric invariant is due 
to the existence of a conical singularity. 

The last two-dimensional cut, namely the h — 7 plane 
at fixed 4>i is significantly more complicated and we have 
not been able to find any simple shape to represent this 
part of the phase diagram. However, using the technol- 
ogy that we develop below for the more easily visualized 
surfaces, we analyze the h — 7 plane in Sec. |III C 
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FIG. 3: (insets) Equivalent graphical representation of the 
phase diagram of the XY model in the 7 — <j> plane, where 
7 G [0,oo) and 4> £ [0,27r]. The right inset shows the para- 
magnetic disordered phase and the left inset represents the 
ferromagnetic phase. It is clear that in the latter case there 
is a conical singularity developing at 7 = which represents 
the anisotropic phase transition. The plots show bulk Eulcr 



integral Xbuik (e) as defined in Eq. 35 demonstrating the jump 
in Xbuik at the phase transition between the paramagnet and 
ferromagnet in the limit e — > + . 



II. UNIVERSALITY OF THE EULER 
INTEGRALS 



We now wish to show that the Euler integrals charac- 
terizing various phases of the XY model are universal to 
such phase transitions due to critical scaling of the met- 
ric. We begin by considering the transverse-field Ising 
(TFI) model with 7 = 1, h € (-00,00), and <p £ [0,2tt). 
For this model, it is knowrP^that the metric tensor, and 
thus the associated curvatures, obey certain scaling laws 
near the QCP. Therefore, since the boundary of the phase 
is at such a QCP, critical scaling theory is encoded in the 
boundary Euler integral. 

However, knowing the boundary Euler integral is suffi- 
cient to determine the bulk integral. To see this, consider 
the region he (— 1 + e, 1 — e) for small positive e. Since 
the region only spans a single phase, there are no ground 
state degeneracies within this region, meaning the h — <j) 
surface is homeomorphic to an open cylinder. Because 
an open cylinder has Euler characteristic 0, the Gauss- 
Bonnet theorem becomes 



Xbuik ^ Xboundary ■ 

boundaries 



(36) 



We want to solve for the bulk Euler integral, in the limit 
that the boundaries of the region are taken to th e p hase 
boundary (e — > + ). However, according to Eq. 36 the 
bulk Euler integral is just minus the sum of the boundary 
integral, which are much easier to solve for. This "bulk- 
boundary correspondence" is what allows us to use crit- 
ical scaling theory to determine the bulk Euler integral 
for each phase. 
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A. Example: Exact metric of the XY chain 

As an initial demonstration of this method, consider 
the exact expression for the metric of the TFI model, 
given in Eq. [31] For a diagonal metric along a curve of 
constant h, the geodesic curvature reduces to 



Ji const. 



y/9<j><f>9hh 



For the case \h\ > 1, this gives 
-4//i 2 



k g {h) = 



= -4(h 2 - I) 1 ' 2 



(37) 



(38) 



/t 4 2 -i) 



Integrating over one of the critical lines, h = ±1 and 
4> e [0, 2tt), gives Xboundary = 1) = 0. To get some 
intuition as to what the boundary Euler integral of zero 
means, consider the three-dimensional embedding shown 
in Fig. [2] A curve with k g = is, by definition, a 
geodesic. This makes sense, since the circle at ft. = 1 
is clearly a geodesic of both the cylinder and the hemi- 
sphere. In general, a smooth curve on a cylindrically- 
symmetric surface will be a geodesic if the radius is at a 
local extremum, i.e., dr/dz = 0. This is clearly satisfied 
in the case of the TFI model, because dr/dh is finite near 
the QCP, while dz/dh^ oo (see Fig. [2j inset). 

Similarly, for the limiting point at h = oo, we can 
calculate the boundary Euler integral: 



Xboundary (h) 



2tt 



k„dl = — 
9 2vr 



2tt 



h 2 - 1 



h 2 



-1. 



k g{ h )\J 94>4>( h ) # 
(39) 



By the Same logic, Xboundary 

Therefore, using Eq 



3(i 



-1 in the limit h - 
we quickly obtain the 1 



0-1 



breakup of the bulk Euler integral for the TFI model. 

To further illustrate the analytical power of this 
method, we can now compute the bulk Euler integral 
of the 7 — <f> plane for arbitrary \h\ < 1, as defined in 
Eq. [35} By a similar analysis as before, one finds that 



Xboundary (7) 



k g{l)\j 9^{l) d<t> 



1 

1 7^0 1 
V 2 (7 + l) ~^ 



(40) 



Since the limit 7 — > 00 corresponds to a geodesic (see 
Fig. [3L giving Xboundary (7 00 ) = 0, we see that the 
bulk Euler integral is xtmik 



1 

V2- 



B. Universality from critical scaling of the metric 

We now use critical scaling theory to find the Euler in- 
tegrals of these phase boundaries for more general mod- 
els. Consider first the case of an arbitrary model in the 



TFI (a.k.a. 2D Ising) universality class. We know from 
the scaling theory of Ref. [19] that the metric diverges 
near the QCP with a power law set by scaling dimen- 
sion of the perturbing operators, d\H . For example, in 
the transverse-field direction, the metric must scale as 
9hh ~ \h-hc\~ 1 for arbitrary models in the Ising univer- 
sality class. Similarly, since the parameter <fi is marginal 
near the Ising critical point, the singular part of g^ has 
scaling dimension +1. Adding in the regular part of g^ 
to get a non-zero value near the phase transition, we see 
that to leading order ~ A + B\h — h c \. Plugging 
this into the formula for the boundary Euler integral, 
one finds that 



const. 



Xbou 



\h - h. 



1 1/2 h- 



0. (41) 



Therefore, the boundary (and thus bulk) Euler integral 
of the Ising phase transition is protected by the critical 
scaling properties of the metric tensor. In terms of the 
geometry of the three-dimensional embedding, adding ir- 
relevant perturbations to the Hamiltonian will shift the 
critical point and deform the shape away from the critical 
point. However, the phase boundary between the ferro- 
magnet and paramagnet will remain a geodesic (dr/dz 
will remain zero). The fact that the geodesic curvature 
is zero on the phase boundary has an intuitive physical 
interpretation. Consider two points on the phase bound- 
ary. The geodesic defines the line of the shortest distance 
between these two points in the Riemannian manifold 
defined by the metric tensor g. It is clear that this line 
should be entirely confined to the phase boundary since 
any deviations from it result in moving toward the di- 
rection of the relevant coupling, along which the metric 
tensor diverges. Since the phase boundary coincides with 
the geodesic, the geodesic curvature is zero by definition. 

To understand the more complicated anisotropic direc- 
tion, we expand the Hamiltonian around 7 = 0. Close 
to the QCP (I7I <C 1), the spectrum is gapless at a sin- 
gle momentum fco = cos _1 (/i) € (0,7r), around which we 
can linearize the equations. Then the linearized mode 
Hamiltonian is 



H k = - 



(k — fco) sin /c 7e i ^sinfc 
7e _l *sinfco —(fc— fco) sin fco 
;o) sin fco 



(42) 



((fc — fco) sin fco) o z k — (7 sin fco) [u^ cos <p — <r\ sin <j>} 



where cr^ x ' v ' z ^ are pseudo-spin Pauli matrices. The pres- 
ence of sin(fco) in both terms suggests fine-tuning, but 
this turns out to be unnecessary. Therefore, we consider 
a more general Hamiltonian of this form: 



Hk' = f}tfa%, + af3-f {<r%, cos + erf, sin < 



(43) 



where a and (3 are arbitrary constants and fc' = fc — fco 
is the momentum difference from the gapless point. This 
linearized Hamiltonian has 



(44) 



The scaling limits of g 77 and are now relatively 
straightforward to compute. The formulas are, as before, 



577 ~ AL Y, 



d9 k > 



9<p<p = ±l^2 sil,2 ( e k') ■ 
k' 

(45) 

In the thermodynamic limit, we convert the sum to an 
integral and define the scaling variable 



K = k' '/7 , 



(46) 



which goes from k = — oo to oo in the scaling limit, I7I <C 
1. Thus, 



1 

8^ 



,dk' 



a 2 7 2 



577 



8tt 

8tt j k 
8717 



-dn 



k , a 2 "/ 2 + k' 2 
017 



dk' 



97 



dk' 



zdn 



1 67 



(47) 



Finally, we use our earlier equation for the boundary Eu- 
ler integral to arrive at 



Xboundary 



472 7 



f 



-r 



-1/2 



V2 



(48) 



Thus, for all models whose low-energy Hamiltonians are 
described by Eq. |43j the bulk Eulcr integral between the 
anisotropic QCP and the geodesic at 7 = 00 remains . 



C. Robustness against angular distortions 

The previous section demonstrated robustness of the 
Eulcr integrals at phase transitions for the case where the 
metric is diagonal. However, our physical intuition from 
the theory of continuous phase transitions suggests that 
this robustness should be even more general, allowing ar- 
bitrary perturbations to the model as long as they do not 
change the scaling properties of the critical point (with 
regards to traditional observables) . Therefore, in this sec- 
tion we demonstrate that perturbations which satisfy this 
constraint while introducing off-diagonal components to 
the metric nevertheless do not change the value of the 
Euler integrals. 

In the h — <fi plane, a simple method for introducing 
off-diagonal terms to the metric is to allow 7 to vary in 
the vicinity of the QCP. Let 7 be some arbitrary function 
"l(h, 4>), with the restriction that 7 > so that we remain 
in the same phase. With this additional freedom, we get 
a new metric g'(h, <j)) such that 

ds 2 = ghhdh 2 + g^d<j) 2 + g^dj 2 + 2g hj dhd r y 

= g'hhdh 2 + g'^dtf + 2g' u dhdct> . (49) 



Noting that dj = (dh^f)dh + (d^^dcj), we find 
9hh = 9hh + 2g hl (d h -i) + g.yjidhj) 2 

g'4,4, = gw + 577(^7) 2 
g'h<ti — g~t~< 



(50) 



Close to the critical point, only one term diverges: g^h ~ 
\h-hc\- 1 => g' hh - \h-h c \- x , while both ^ and 5 ^ 
remain finite near the critical point. Thus, g' is asymp- 
totically diagonal near the critical point, our earlier ar- 
guments still work, and the boundary Euler integral re- 
mains zero. 

Not surprisingly, the non-integer bulk Euler integral of 
the anisotropic phase transition is more sensitive to de- 
tails of the perturbation. For instance, the most naive 
option of giving the transverse field a functional depen- 
dence (h —> h(-f, <f)) changes the value of the bulk Euler 
integral, which is not surprising given that h is a relevant 
perturbation at this phase transition. This can be traced 
back to the fact that modifying h changes the position of 
the gapless momentum ko (see Eq. 42 ) , strongly affecting 
the low-energy physics near the critical point. 

In the absence of physical parameters to modify, we 
instead consider modifications to the low-energy Hamil- 
tonian. In particular, consider a slightly more general 
Hamiltonian of the form: 

H k > = j3k' al, + (3j [a x (4>)a k , cos0 + a y (4>)cr y k , sin0] , 

.(51) 

where we demand that the functions a x ,y are periodic 
(a XjV (0) = a X:V {2-K)) and positive, such that the az- 
imuthal Bloch angle still wraps the sphere once as we 
take (f> from to 2ir. 

To determine if the Euler integral Xbomidary = — 1/ v2 
is protected, we numerically solve for the boundary Euler 
integral for a variety of functions a x ^ y . In doing so, we re- 
quire an additional constraint to ensure that this integral 
is well-defined: the metric must be positive definite, i.e., 
its determinant g — EG — F 2 must be non-zero. We have 
tested a number of functions satisfying these constraints, 
and found that all of them have Xboundary = — 1/V2 as 
expected. Given that the most complex functions we 



tested (a x = 1 



and a y 



have 



no special symmetries, we postulate that the Euler in- 
tegral is identical for all functions satisfying the above 
constraints; however, we are unable to analytically prove 
such a statement at this time. 



III. CLASSIFICATION OF SINGULARITIES 

Using scaling arguments, we have demonstrated the 
robustness of the geodesic curvature and the bulk Euler 
integral for situations where the boundary of the param- 
eter manifold coincides with the phase boundary. One 
obvious difference between the model in the transverse 
field {h — 4>) and the anisotropy (7 — <j)) planes is integer 
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vs. non-integer values of the Euler integrals. In this sec- 
tion we show how this difference comes from the nature of 
the singularities at the respective phase boundaries. We 
identify two types of geometric singularity: integrable 
singularities, as in the case of the h— <f> plane, and coni- 
cal singularities, as in the case of the 7 — <\> plane. Finally, 
in the h — 7 plane, we identify a third type of singularity, 
known as a curvature singularity. We discuss general con- 
ditions under which these singularities should occur and, 
for the case of conical singularities, identify the relevant 
parameters in determining the boundary Euler integral. 



A. Integrable singularities 

A simple question which we must ask before classify- 
ing the geometric singularities of the XY chain are what, 
precisely, do we mean by singularities? A simple defini- 
tion, namely the divergence of one or more components 
of the metric tensor, is certainly a useful tool for diag- 
nosing the presence of phase transitions in practice^l2U. 
However, we claim that this singularity is less fundamen- 
tal from a geometric standpoint. For instance, in the 
case of the TFI model, the transverse field component of 
the metric tensor diverges as g^h ~ \h — 1| near h = 1. 
However, this divergence can be removed by simply repa- 
rameterizing in terms of h! = v h — lsgn(/i— 1), for which 
9h'h' ~ 1- Therefore, we need to look elsewhere for infor- 
mation about the fundamental nature of the singularities 
in the quantum geometry. 

Since the issue with the metric tensor was its coordi- 
nate dependence, natural quantities to look at are the 
geometric invariants introduced in Sec. |I A[ which are 
coordinate-independent. Of these, the Gaussian curva- 
ture is the obvious choice^. We therefore classify sin- 
gularities here and in the rest of the paper based on the 
Gaussian curvature, K, and its invariant integral, xtmik- 

For the case of the TFI model, the curvature does not 
diverge near the critical point. This can be easily seen in 
the equivalent three-dimensional manifold (Fig. [2]), where 
the curvature goes from that of a cylinder (K = Qj to that 
of a sphere (K = 1/a 2 , where a is the radius), both of 



which are finite. As we show more explicitly in Sec. Ill C 



one can derive this non-divergent result by using the scal- 
ing forms of the metric tensor to get K ~ const. 

However, critical scaling theory does not demand that 
the curvature is a smooth function of the transverse field. 
Indeed, we expect it to be singular (like most other quan- 
tities) in the vicinity of a phase transition, which mani- 
fests in the TFI chain as a jump of K between the fer- 
romagnet and the paramagnet. However, the curvature 
is finite at all points, and is therefore completely inte- 
grable when determining Xbuik- Therefore, we refer to 
these jumps in the curvature as "integrable" singulari- 
ties. We note that visually these integrable singularities 
correspond to points where the manifold changes shape 
locally, but in such a way that the tangent plane evolves 
continuously, so that no cusps or other points of curva- 



ture accumulation occur. 



B. Conical singularities 

The anisotropic pha se tran sition at 7 = is an example 
of a conical singularit} * 24 * 31 !, which can easily be seen in 
Fig. [3] While the specific value of 1/ \/2 for the bulk Euler 
integral is likely specific to this particular class of models, 
we claim that the existence of conical singularities is in 
fact a much more general phenomenon. 

More specifically, we expect conical singularities to oc- 
cur in situations with two inequivalent directions orthog- 
onal to a line (or a higher dimensional manifold) of crit- 
ical points, as long as the orthogonal directions have the 
same scaling dimension.^ Denote these directions Ai and 
A2, with the critical point at Ai = A2 = 0. At 7 = in 
the anisotropic XY model, the parameter <p has no ef- 
fect on H, so in the FM phase this model satisfies the 
criteria for a conical singularity with Ai = 7 cos 4> and 
A2 = 7 sin <j). 

For simplicity we also assume that the metric has cylin- 
drical symmetry, as in the case of the anisotropic tran- 
sition in the XY model. In the previous section we ver- 
ified numerically that this singularity, and in particular 
the boundary contribution to the Euler characteristic, is 
protected against breaking of the cylindrical symmetry. 
We nevertheless use this assumption to simplify our anal- 
ysis. To ensure cylindrical symmetry, the metric tensor 
should be diagonal in the A — <j) plane with the leading or- 
der asymptotic of the diagonal components of the metric 
tensor scaling as some power laws: 



9\\ 



AX 



9U 



(52) 



where A and B are arbitrary positive constants and we 
generally expect a > 0, f3 > 0. However, if we define Ai — 
A cos 4> and A2 = A sin <£>, then the demands of uniform 
scaling place an additional constraint on the values of the 
scaling dimensions. To see, this consider the components 
of the metric in "Cartesian coordinates" : 



9u 



9x\ 



<9A7 



■ 9<t><t> 



9\\ 



A 2 



9<j>4> 



Af 
A 4 



AX~ a cos 2 . 



522 
gi2 



= A\- 



- B\?- 2 sin 2 , 
BX P-2 „„„2 



' cos 



= (AA^ + BA^^cos^sin^ 



(53) 



We clearly see that the scaling dimensions of Ai and A2 
are the same at all angles <p if an d only if the exponents 
satisfy the relation: 



p = 2 



(54) 



Note that the condition j3 > now means that < a < 2. 
The constants A and B are non-universal, but we ex- 
pect that their ratio B/A - which defines the anisotropy 
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FIG. 4: Graphical representation of the metric in the Ai- 
A2 plane, where a conical singularity (critical point) is at 
Ai = A2 = 0. Panel (a) shows the simple case of the XY 
model in a transverse field, for which a = ft = 1, B/A = 2, 
and the metric is cylindrically symmetric about the origin. 
Panel (b) shows the more general case with Hamiltonian given 
by Eq. |51[ in which the cylindrical symmetry has been bro- 
ken by using the functions a x (cj>) = 1 + ~ cos(2</>) — | sin(2</>) 
and a y (({>) = 2 + sin(4>). While the cylindrical symmetry is 
gone, we numerically find that this model has the same con- 
ical angle as when a x = a y — const, yielding the same bulk 
Euler integral. The metric is plotted in both panels by show- 
ing the "shape of the circle" near each point, i.e., the blue 
ellipses show contours of constant radius in the A1-A2 plane. 
Since the metric diverges near the critical point, the size of 
the ellipses gets smaller. 



of the metric tensor - will be a universal number for a 
given class of models. For the anisotropic transition this 
ratio is B/A = 2 and the exponents are a = = 1 (see 
Eq. 47). Interestingly the point h — 00 in the h — cj> plane 
(corresponding the spherical cap - see Fig. [2]) also has 
the form of a conical singularity if we use Ai = \ cos(0) 
and A 2 = \ sin(0) with a = 0, = 2, and B/A = 1. 




These exponents describe a non-singular point in param- 
eter space with cylindrical symmetry. 

Given a conical singularity, we can now easily find the 
Euler integral. Using the same formulas as earlier for the 
case of cylindrical symmetry, 



Xboundary , (55) 

V.9AA 

/rv \ f~B 

(56) 

Using this formula for the anisotropic phase transition 
of the XY model yields Xboundary = —l/y/2, as found 
earlier. For a demonstration of the contours of the metric 
in this model, see Fig. |1J 

Let us point that if we have a non-singular point like 
h — > 00 at fixed 7, for which a = and f3 = 2, we have 
an additional requirement that gn = (722 and 312 = 0, 
implying that A = B. This follows from the fact that the 
metric must remain regular at A = 0. Then from Eq. [56] 
we find that the boundary contribution of the isotropic 
point will be Xboundary = — 1. In a three-dimensional 
embedding, this indeed looks like a hemisphere, which is 
non-singular. Therefore, for such smooth "singularities," 
the manifold is also guaranteed to be locally equivalent 




FIG. 5: Three dimensional plot of K(h,y), similar to that 
found in Ref. 1111 The graph shows curvature singularities 
along the two phase transitions h = 1 (only positive h are 
shown) and h < 1, 7 = 0. Unlike the Ising phase transition, 
the curvature singularities near the anisotropic phase tran- 
sition are non-integrable, leading to divergent Euler integral 
within each phase. 



to the hemisphere. 

Finally, for the case of a multi-critical point, one can 
try to apply similar logic. However, due to the asym- 
metry in the scaling dimensions, along some direction 
the metric will be infinitely anisotropic near the critical 
point. This infinite anisotropy is not generally remov- 
able by rescaling the couplings. Therefore, the conical 
singularity breaks down and the curvature can become 
non-integrably singular. 



C. Curvature singularities 

If we now consider the third two-dimensional cut of the 
XY model, namely the h — 7 plane, we find that the cur- 
vature displays a number of additional singularities. The 
structure of these singularities can be seen in Fig. [5j It is 
clear from the plot that as expected there are singulari- 
ties near the two phase transitions: integrable singulari- 
ties near the Ising transition (\h\ — 1) and non-integrable 
singularities near the anisotropic transition (\h\ < 1 and 
7 = 0). These singularities meet near the multi-critical 
point (| ft. I = 1 and 7 = 0) resulting in a very singular and 
non-monotonic behavior of the curvature. 

While, unlike the h — and 7 — <j) planes, there are 
no obvious finite protected Euler integrals in the h — 7 
plane, the exponent of the curvature divergence can be 
found from scaling arguments. For instance, near the 
Ising phase transition {\h\ = L 7 > 0), the simple argu- 
ments of Venuti and Zanardf^ indicate that the metric 
components diverge as guh ~ \h — ft c | _1 , 577 ~ 1, and 
gh-y ~ 1. The non-divergence of g 77 is due to the fact 
that 7 is a marginal parameter, such that the scaling di- 
mension of the singular part of <? 77 is +1. However, there 
is also a non-zero non-singular part, which is the leading 
order term near the critical point. Similarly, the scaling 
dimension of gh-y is zero, such that there is a jump singu- 
larity near the critical point. Plugging these divergences 



into Eq. 11 and assuming a smooth dependence on 7 as 
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long as we are far from the multi-critical point, we find 
that K ~ 1, which matches with the jump singularity in 
K found at the Ising phase transition. By contrast, near 
the anisotropic phase transition (7 = 0, \h\ < 1), both 
the transverse field h and the anisotropy 7 are relevant 
parameters, with scaling dimension —1. Therefore, the 
metric components scale as g^h ~ 9h-y ~ 577 ~ \l\ ■ 
By the same logic as before, this leads to a divergent 
curvature with K ~ I7I ■ 

Finally, we point out some subtleties regarding the cur- 
vature far from critical points. First, while the curvature 
does not diverge near h = 00 or 7 = 00, it is not immedi- 
ately apparent whether the bulk Euler integral diverges 
in this limit. Therefore, in Appendix |Bj we show that 
a simple reparameterization allows one to make both K, 
yjg, and the limits of integration simultaneously finite, 
except near quantum critical points at finite h and 7. 
Therefore, it becomes clear that the divergence in the 
Euler integrals comes strictly from the curvature singu- 
larities near the quantum phase transitions. 

Second, we note that the metric component gun van- 
ishes at 7 = for all h > 1, since the ground state is the 
fully-polarized spin-up state along this entire line. The 
general intuition from Riemann geometry is that if the 
determinant of the metric vanishes, then the curvature 
diverges, since the determinant appears in the denomi- 
nator of the curvature formula (Eq. 11 ). However, as we 
show in the next section, the curvature does not in fact 
diverge for fundamental physical reasons. 



force M., which is generally non-degenerate. By the argu- 
ment above, along the symmetry line the ground state of 
H will be independent of A until the gap in the Hamilto- 
nian closes, e.g., until the system undergoes a quantum 
phase transition. Thus the metric near this symmetry 
line will be singular with a vanishing determinant.^ 

We now investigate why, despite this "singular" metric, 
the curvature nevertheless remains analytic in the vicin- 
ity of such a line of symmetry. Assuming we are far from 
any critical points (i.e., with a gapped spectrum), the 
components of the metric tensor near the fully polarized 
state should be analytic and can be written as 



S 2 fxx 
Sfsx 



Sfsx 

j 66 



(58) 



All components /y are smooth functions of A and 8. For 
small 8 they can be approximated as being independent 
of 8, fij — fij{h), since the leading asymptotic of the 
curvature in the limit 8 — > will be determined by the 



explicit dependence given by Eq. 58 



Using the explicit expression for the curvature (Eq. 11 ) 
and counting powers of 8, we see that the only possible 
divergent term in the curvature as 8 — > is given by: 



K 



1 



d 



6*VJd\ V/aa 



X6 



(59) 



where / = det(/y •) and 



D. Metric singularity near lines of symmetry 

Consider the line \h\ > 1, 7 = in the XY model. 
The ground state along this entire line is fully polarized 
along the direction of the transverse field, which is clearly 
the ground state at h = 00. Then, since along this spe- 
cial XY-symmetric line the Hamiltonian commutes with 
'-'total' the fully polarized eigenstate must remain the 
ground state until a gap closes. Note that this argument 
continues to hold even in the presence of integrability- 
breaking perturbations, as long as the z-magnetization 
remains a good quantum number. Therefore, such a line 
of unchanging ground states and thus vanishing metric 
determinant is a robust feature of this class of models. 

More generally, one can create such fully polarized 
ground states by considering a family of the Hamilto- 
nians 



H(X,8)=H (X,S)-XM , 



(57) 



where 5 is a symmetry breaking field such that Ho and 
M. commute at 5 = for any value of A. Physically M. 
represents the generalized force (z-magnetization in the 
above example) which is conserved at the symmetry line. 
In general M. can also depend on A and 8 as long as Ho 
and M. commute at 8 = 0. Clearly in the limit A — > 00 
the ground state of the Hamiltonian is the fully polarized 
state (the state with largest eigenvalue) of the generalized 



1 xs - 



1 



2.9 

Therefore 



ffAA" 



dgss 
dX 



9sx- 



dg. 



XX 



08 



fxx 
2/ 



dfss 
dX 



-Vsx 



K 



d 



2<5 2 V7<9A 



dxfss — 

77 



(60) 



(61) 



In order to see that the curvature is not divergent near 
this line of symmetry, we need to show that the following 
difference vanishes 



dxfss - 2fsx = dxgss - Zdsgxs = 0. 



(62) 



This is indeed the case for such general lines of sym- 
metry, as we prove in detail in Appendix [X] Therefore, 
the singular term in the curvature vanishes, and thus the 
curvature does not diverge near the symmetric line. 

A more geometric view of the absence of a curvature 
singularity for metrics satisfying Eq. [62] can be seen by 
mapping this metric to the surface of the cone in either 
Euclidean or Minkowski space. By matching the metric 
in Eq. [58l 



ds 2 = S 2 fxx dX 2 + 8 dxfss dXdS + f u dS 2 , 



(63) 



to that of a surface in cylindrical coordinates, 

ds 2 = dr 2 +r 2 dc/) 2 + dz 2 , (64) 
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and looking for the solution of the form 

r = 6, z — 5z(X), <t> — (f)(\) 

we find 



66 



1. 



r/A 



rfA 



(65) 



(66) 



Clearly this solution can only work if fgg > 1. In the 
opposite case of fss < 1 we can match the metric to a 
surface in Minkowski space 



ds 2 = dr 2 + r 2 d<b 2 - dz 2 . 



with 



dz N 2 



(67) 



(68) 



In both cases the equivalent surfaces lie on cones z oc 
r oc S, and thus have vanishing curvature. There is fur- 
ther freedom of reparameterization in the conical map- 
ping (Eq. 65) by, for example, choosing z — CSz(X) 
and r = CS, with C an arbitrary positive scaling fac- 
tor. This reparameterization changes the condition of 
switching from the Euclidean to the Minkowski mapping 
from fss = 1 to fgg = C 2 . However, in all the cases 
the mapping remains locally conical because z oc r, and 
thus gives vanishing curvature. Finally, we note that the 
curvature generally becomes finite - and analytic - once 
one takes into consideration non-singular sub-leading (in 
S) terms in Eq. 58 



IV. MEASURING THE METRIC 

Formally the metric tensor can be expressed as a stan- 
dard response function^ and thus is in principle measur- 
able. However, unlike the Berry curvature, which natu- 
rally appears in the off-diagonal Kubo-type response^, 
the metric tensor a ppear s either as a response in imag- 
inary time dynamic d 16 ! 17 ! or as a response in dissipative 
systems using specific - and usually not physically jus- 
tified - requirements for the dissipation^. However, 
as shown in Ref. [TS] for the specific situation of non- 
interacting particles, the geometric tensor characterizing 
the Bloch bands can be measured through the spectral 
function of the current noise. Here we extend this idea 
to arbitrary systems and couplings. 

The geometric tensor can be represented as (cf. Eq. 
(B8) in Ref. [16| 



2ir oj 2 



(69) 



where 



dte~ i " t {d ll H{t)d v H^)) c (70) 



is the Fourier transform of the connected ground state 
non-equal time correlation function of the generalized 
forces, a.k.a. the noise spectral density, and 9„iJ(t) = 
e lHt d ll H(Q)e~ lHt is the generalized force in the Heisen- 
berg picture. The metric tensor is the symmetric (real) 
part of the geometric tensor and thus can be expressed 
through the symmetrized spectral density. Following the 
insight of Neupert et alP^ we interpret S^(lo) as the 
Fourier transform of the non-equal time noise correlation 
function of two generalized forces, which is relevant ex- 
perimentalljJ231. For example, in mesoscopic systems, the 
current noise spectrum S^ v {oj) can be measured in shot 
noise experiments, where current J M corresponds to the 

suggests a simple and 



(in 



generalized force Qa^H- Eq. 
general way of measuring the metric tensor in interact- 
ing many-body systems by analyzing equilibrium noise. 
We believe that, for sufficiently large systems, such sym- 
metrized noise correlations should be measurable with 
negligible effects of measurement back-action. 

While the method of measuring the metric tensor 
through noise is conceptually simple, it cannot be eas- 
ily implemented in systems such as cold atoms, where 
measurements are often destructive. Below, we discuss 
two real-time protocols which offer the possibility of ob- 
serving the metric tensor via destructive (single-time) 
measurements. We then mention an additional protocol 
involving instantaneous quenche s of th e external param- 
eters. We note that both ramp^MH and quenches^^U 
are routinely achieved in isolated cold atom systems. 

Consider performing real-time ramps of some parame- 
ter A M in a gapped system, starting from the ground state 
at the starting point \. It has been shown elsewhere^ 
that, for a square root ramp with A^(i) — A^(i/) = 

1/2 

[v(tf — t)] , the leading order correction to the energy 
in the limit v — > is given by 



(H) =E + vg^ + O(v 2 ) 



(71) 



where H is the Hamiltonian, Eq is its ground state en- 
ergy, and g^ is diagonal component of metric along the 
ramping direction. 

However, the square root ramp is singular near tt, and 
therefore may be difficult to implement. We now show 
that the metric can also be measured via a more eas- 
ily implemented linear ramp, at the cost of requiring a 
harder measurement: the quantum energy fluctuations. 

Consider a linear ramp A„(t) — X^(tf) — v(t — tf). 
From Ref. 1391 we know that the wave function at tf will 
be given in its instantaneous eigenbasis by 

\il>) = (1 + /3w 2 )|0) -iv^2 a n \n) + 0(v 2 ) , (72) 

n#0 



where a r 



jS? and ^ = -(V2) Kl 2 * which 

serves to keep the wave function normalized up to or- 
der v 2 . The energy fluctuations are given by AE 2 = 
(4>\H 2 \4>) — (ip\H\tp) 2 . Without loss of generality, we may 
offset the Hamiltonian such that the ground state energy 
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is Eq = 0. Then, up to order v 2 



AE 



v 2 K 
v 2 \a r 

„.2 \ " 



(n\H 2 \n) 



\ 2 E 2 



i\d^H\Q){Q\d^H\ 



E 2 



(73) 



Therefore, by measuring the energy fluctuations for 
different ramp rates and extracting the leading order 
(quadratic) term, we can extract diagonal terms of the 
metric along a given direction. Let us point that if we 
start the ramp in the ground state then the energy fluc- 
tuations are equal to the work fluctuations, so the metric 
tensor can be extracted by measuring work fluctuations 
as a function of the ramp rate. 

A third possibility for measuring the metric tensor is 
by measuring the probability of doing non-zero work for 
small quenches in parameter space. This is in some sense 
true by definition: if \ipo(X)} is the ground state manifold, 
then the probability of doing zero work (i.e. ending up 
in the ground state) after a quench from A M to A M + dX^ 
is just 



P(W = 0) = \(M^)\M^ + dX^))\ 

= 1 - g^dXl ; 
P(W ? 0) = g^dXl . 



(74) 



As noted elsewhere, this quantity is equivalent to the 
time-averaged return amplitude G(tJ^t 



P(W = 0) 
G{t) 



km i / G(t)dt , 

{MK)W H{ ^ ]t e- lH(K+dK ' ]t \M\d) 



which is related to the well-known Loschmidt echo L(t) 
by 



Ht) = \G(tf 



(75) 



The Loschmidt echo is also the probability of returning to 
the ground state (doing zero work) after a double quench 
of duration t from A p to A M + dX^ and baclPSl. While 
energy distributions and the related Loschmidt echo are 
in principle measurable by a variety of methods, we note 
that there has been important recent progress in propos- 
ing measureme nts of these quantities using few-level sys- 
tems as a probe^L^I. 

Finally, we point out that one can reconstruct the 
full metric tensor solely from measurements of its diag- 
onal components. Consider a two-parameter manifold 
{X Xl X y ). First, measure the diagonal components g xx 
and g yv using one of the procedures described above. Sec- 
ond, measure a specific off-diagonal element by varying 
X x and Aj, simultaneously. For example, if we define the 
variable A^ = (X x + X y )/2 and ramp or quench along the 



line X x = X y , we can obtain g ww . Finally, noting that for 
this protocol dX x = dX y = dX w , we see that 

ds 2 = g ww dX 2 w = g xx dX 2 x + 2g xy dX x dX y + g yy dX 2 
— gwwdX 2 = (g xx + 2g xy + g yy ) dX 2 

gww gxx gyy 



9xy 



(76) 



This procedure can be easily generalized to an N- 
parameter manifold by performing pairwise measure- 
ments using a similar tricks as above. 



V. CONCLUSIONS 

In conclusion, using the quantum XY model as an ex- 
ample, we have analyzed the Riemann manifold of a sim- 
ple ground state phase diagram. We identified a new ge- 
ometric characteristic - the bulk Euler integral - which 
characterizes the phases of matter. Based on the value of 
this Euler integral, either integer, non-integer, or unde- 
fined, we have classified three types of singularities in the 
Gaussian curvature: integrable, conical, and curvature 
singularities. We showed that integrable singularities oc- 
cur for phase transitions where one parameter is marginal 
or irrelevant while the other is relevant. Similarly, conical 
singularities emerge when the phase transition occurs at a 
single critical point with two "orthogonal" relevant direc- 
tions that have the same scaling dimensions. And finally 
near the multi-critical point with two inequivalent rele- 
vant directions we found curvature singularities which, 
similar to black holes, are non-removable non-integrable 
singularities in the quantum metric space. Finally, by in- 
troducing additional techniques for measuring the metric 
experimentally, we point out that this geometric informa- 
tion should be experimentally accessible. 
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Appendix A: Proof of Eq. [62] 



We now prove that Eq. 62 indeed holds for the metric 
given by Eq. [58] near the symmetric line. We will rely 
on the fact that, at this symmetric line, the ground state 
does not depend on A, so 9a |0) = 0. Thus, 



14 



dxgss = d x ]T(0|^|n)H%|0) = £ (<0|§>> + <ofe|rc}) (n|S 5 |0) + (0|£|n) (<nfo|0) + (n|^|0) 



(Al) 



n#0 L 



E 



<0|S|n}<n|S A |0} + <0|^|n}<n|%|0) 

(<0|Sf^|n) + (ofo|n)) (n\d s \0) + (0\K\n) ((nfe|0) + (n|0^|O) 

Therefore, we now find that Eq. [62] holds: 

dxg 5S -2d s9sx = Y,[( &dx\n)(n\ds\0) + (0\^s\n)(n\^ds\0) 

[{0\t s \m){m\dx\n){n\d s \0) + {0\fr s \n){nfc\m){m\d s \0) 



(A2) 



. 



(A3) 



The last equality follows by observing that if we inter- 
change indices n and m in the second term in the last 
sum, we get a term that exactly cancels the first one. 
This follows from 



\K\n) 



(m\d\\n). 



(A4) 



In addition, the m = term vanishes because 

(0\d x \n) = -(0\K\n) = . (A5) 

Appendix B: Choice of parameters 

The goal of this section is to show by example that, if 
the bulk Euler integral 27TXbuik = J KdS diverges, then 
that implies that the curvature must diverge at some 
point. This is not a priori obvious, because the invariant 
area dS — y/gd\id\2 can also diverge, either because the 
metric diverges or because the metric is finite but the 
parameters Xi have infinite range. We show that, for the 
h — 7 plane of the XY model, such divergences can be 
removed by a suitable choice of coordinates. 

One natural thing to attempt to do is the go to "unit- 
less" coordinate systems, dh — > Jg~hh,dh and d^y — > 
y/g-y 7 dj, in which a one-parameter metric g\\ would be- 
come flat (i.e., become A-independent). While this does 
not quite work the same for two parameters, since ghh 
depends on 7, we nevertheless use a variant of it below 
to get a more well-behaved metric. As we will see, the 
new parameters have a finite range, and the metric is 
much more well-behaved. 

Consider first the case of the transverse field h. We 
wish to define a new parameter £ such that d£ = ^Jghhdh. 
The leading order h dependence of ghh is ghh ~ (1-h 2 ) m 





FIG. 6: Geometric invariants of the XY model in the £-77 
plane. The three dimensional plots show Gaussian curvature 
K (a) and metric determinant g (b) for the entire (finite) 
range of parameters £ and 77. As discussed in the text, the 
invariants clearly only diverge near the critical line at 7 = 0, 
except for the determinant g, which has an integrable diver- 
gence near £ = 77 = n. 



the ferromagnet and ghh ~ h 2 (h 2 -i) m ^ ne P arama g ne f ■ 
Integrating these expressions gives the natural choice 



sin£ 



h if \h\ < 1 
1/h i£\h\>l 



(Bl) 



with the quadrant of £ chosen such that |£| € [0, 7r/2] 
if \h\ < 1 and |£| G [7r/2,7r] if \h\ > 1. We similarly 
reparameterize the 7 direction in terms of 77 such that 
dri = ^ +7) rf7, giving 



(B2) 



77 = 2 tan (V|7|)sgn7 



The range of our new, auxiliary variables is £ G (— 7T, 7t), 

77 e [o,7r). 

Within the FM phase, the metric takes on a simple 
form: 



1 

Off = — cot' 

y« 16 



(I) 



1 

16 



(B3) 
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By inspection, the metric and its determinant clearly 
only diverge at the anisotropic phase transition, which 
is at r\ = in the new parameters. Within the PM 



phase, the metric remains quite complicated. However, 
one can easily calculate the curvature and determinant 
of the metric, which are given by 



K = 8 



csc 2 (£)sec 4 (|)tan 6 (|) 



Scsc(£) 

c 2 (£)+tan 4 (f) 

-1 + 2 esc 2 (£) + tan 4 ( \ ) - 2 esc (£) J- 1 + esc 2 (£) + tan 4 ( 



256 (tan 4 (f) - LJ (-1 + csc 2 (£) + tan 4 (n))' 



(B4) 



These two quantities are plotted for the entire phase di- 
agram in Fig. [6j 

Clearly we have almost achieved our goal, in that the 
invariant area component of the bulk Euler integral, dS = 
y/g d^drj is finite except near a few select points, namely 

• In vicinity of the critical line at ?/ = (i.e. 7 = 0), 
where the curvature also diverges. 

• Near the points \h\, I7I — > 00, which correspond to 
|£| = \q\ = 7r. Here the curvature does not diverge, 
so we expect the divergence of the metric to again 
be removable. 

To see how to remove the divergence in g near r/ = £ — 
ir, we need to understand the asymptotics of g near this 
point. We can do leading order asymptotic expansions of 
the numerator and denominator about this point. Then 
if we define 



1 

csc(C) 



1 

tan 2 (77/2) 



7T — 77 
2 



. (B5) 

we find that the determinant is asymptotically equivalent 
to 



'[2i 



256v~ 



(l) [2(S) 2 + 1-2(S) V 1+ (S)' 



. 2 n 2 



25541 + 0D 

Rewriting this in circular coordinates u 



(B6) 



(B7) 



rcos(9) and 



v = rsm(6), the expression becomes = ^g^/^), 

where the notation g{^ iV } is meant to reiterate that this is 

9vi 



the determinant of the matrix 



The function 



f{6) = sin [sinO - if 



(B8) 



is defined over the interval 9 G [0, 7r/2]. Then the invari- 
ant area is (using the expansions of u and v from above) 



dS 



' sin^l-sin^f-d^f- '''' 



256r 



2^ 




32r 
1 — sin ( 



sin 6(1 — sin ( 



(rdrdO) 



1 



4r sin 9 I 



I dudv 



32 



-drd9 = y/g{ r _o}drd9 . 



(B9) 



So we come to our final result that, by choosing a local 
parameterization (r, 9) as described above for the points 
near |£| = \r}\ = ir, the metric determinant g{ r ,e} is non- 
divergent. We conclude that, after a suitable choice of 
local reparameterization, the invariant area term and its 
integral can be made finite unless curvature diverges. 
Therefore, all divergences in the bulk Euler integral of 
the h — 7 plane are due to the divergent curvature near 
7 = 0. 
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